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Abstract 

The deuteron spectrum in n+p —> 2 H (tt-k) at k n = 1.88 GeV/c and 0$ = 0° 
is explained by considering a AA excitation as the dominant reaction mecha- 
nism for the 2ir production. We present a new theoretical approach based on 
a coupled channel formalism which allows to include the residual interaction 
within the intermediate AA and A./V systems. The corresponding interaction 
potentials Vaa and Van are adopted from a meson exchange model with ir, 
p, u, and a exchange taken into account. The influence of the residual in- 
teraction on the deuteron spectrum is studied. We also predict the angular 
distribution of the two pions. It is shown that this distribution is closely 

connected to the spin structure of the AA excitation. 
13.75.-n, 25.10,+s, 14.20. Gk, 24.10.Eq 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The ABC effect, a mesonic structure of isospin 1 = 0, was first observed by Abashian, 
Booth and Crow in the reaction p + 2 H — > 3 He+(7T7r) [[[]-§]. Later on, it was also shown 
to appear in several other hadronic reactions, like e.g. in n + p —* 2 H +X and in 
2 H+ 2 H— > 4 He +X In all these cases, the ABC structure shows up as a cross section 

peak at missing masses of about 300 - 350 MeV, and with a position and width that varies 
quite rapidly according to the kinematic conditions. The only possible explanation for the 
effect seems to be an enhancement in the double pion production which is caused by a 
strongly energy dependent production amplitude. 

Since the ABC effect is present in reactions of different type, it seems to be quite natural 
to assume that there exists a common underlying mechanism for the 27r production which 
should be able to explain all the different experiments. From a theoretical point of view, 
n + p -^ 2 H (ittt) has to be the reaction to analyze first since it involves only two-nucleon 
states beside the pions and can thus be treated in the most accurate and direct way. For 
this reaction, the ABC enhancement corresponds to a beam momentum of approximately 2 
GeV/c in the laboratory system and to a center of mass energy of ~ 2Ma, i.e. twice 
the mass of the A(1232) resonance. This fact strongly suggests that the two pions are 
dominantly produced via a double A excitation. Indeed, as it has already been shown by 
Bar-Nir et al. [p|,|ll|, the main features of the experimental spectrum can be understood by 
assuming the AA mechanism to be the most efficient one for 2ir production in the energy 
region under consideration. 

In the present work, we will present a new calculation of AA and AiV excitations in 
the n + p — > 2 H (tttt) reaction which explicitly includes the effects of residual interactions in 
the intermediate resonance states fl2yi3|| . Our aim is to show that the 2ir production may 
serve as an interesting tool for further examination of direct AA and AN potentials which 
are not very well known. For this purpose, we will also examine how the spin structure 
of the AA excitation influences the angular distribution of the two pions. Our model is 
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based on a coupled channel approach in a non-relativistic framework. The potentials of 



the residual interactions are adopted from a meson exchange model |H]I5|| . In this work, 
contributions from 7T, p, u>, and a exchange are included. The A resonance is treated thereby 
as a quasi-particle with a given mass and an energy-dependent, intrinsic width. 

Experimental results for the n+p —* 2 H (ittt) cross section in the ABC energy region have 
been obtained by Plouin et al. @] at k n = 1.88 GeV. In order to allow for a direct comparison 
with our theoretical results, the contribution from one pion production has been subtracted 
from the spectra. It has been stressed that there might be some additional background from 
n production [16| due to momentum spread of the beam, but since there is no microscopic 
model for this contribution available, we decided not to perform any further corrections. 
One should keep in mind, however, that not all of the uncertainties have been included in 
the experimental error bars. 

In section [II], we will present the theoretical framework of our model. The results of our 
calculation are presented and discussed in section [IT| The paper concludes with a summary 
in section [TV]. 

II. THEORETICAL FRAMEWORK 
A. Cross section 

We are interested here in the calculation of the deuteron spectrum for the n+p — > 2 rl (tttt) 
reaction, i.e. in the double differential cross section d 2 a / 'dk^dfld- Using relativistic kinemat- 
ics, the cross section is given as 

2E n E n M„ M n d 3 k,i d 3 k 



da = I = ~ET ~ET 77, 77 » /OF \ ( c\ T-i \ S(E P + E n -E d - E nw ) M 1 

I\(s, M p , M n ) E n E p (2tt) 3 (2ir) 3 (2E a )(2E b ) 



where k = \{k a — kb ) denotes the relative momentum of the two pions and E n7T = E a + Eb 
is the sum of their energies. The indices n(p) and d refer to the neutron projectile (proton 
target) and the deuteron ejectile, respectively. The center of mass momentum of the 2ir 
system K = k a + kb is fixed by momentum conservation and equals k n + k p — k^. As usual, 



the function A is defined to be A(s, M 2 , M 2 ) = [s - (M p + M n ) 2 } [s - (M p - M n ) 2 }. Finally, 
| M. | 2 denotes the transition amplitude for the reaction; in the case of the unpolarized cross 
section, an average over the initial and a sum over the final spin orientations has to be 
performed. 

From eq. (|I|) one easily obtains the deuteron spectrum 
d 2 a 1 Mr, .o r k 2 dk 



n #* 7ihr^-M\M\ 2 ( 2 ) 



dk d dQ d (2vr)5 k^ d J AE a E b dE 

Since we are going to calculate the transititon matrixelement in a non-relativistic approach, 
the above expression has to be evaluated in the deuteron rest frame (drf ) where the deuteron 
is well described by the usual non-relativistic wave function ||15|| . Afterwards, the transfor- 
mation of the spectrum into the laboratory frame (lab) is performed by using the simple 
relation 

/ d 2 a \ Ef\k 1 ^) 2 ( d 2 o 



dk d dn d ) lah El? h (k^) 2 \dk d dn a 



(3) 



B. Reaction mechanism and evaluation of the transition amplitude 

In our microscopic model for the n + p -h> 2 H (tttt) reaction, we take into account two 
different reaction mechanisms corresponding to the Feynman diagrams presented in fig. [I]. 
The first mechanism (a) involves the excitation of a double Delta (AA) intermediate state 
where both A resonances subsequently decay into a pion and a nucleon. In the case of 
the second mechanism (b), each of the two pions couples to the same A, hence only one 
resonance (i.e. a AiV system) is excited. 

For the explicit calculation of the two matrix elements, wa take advantage of the principle 
of detailed balance. This allows us to assume that the deuteron is in the initial state and 
absorbes the two pions, which finally will lead to a n + p state. The absorption of only one 
pion with momentum k n is described by the operator 

Fl{K) = e^ 2 (S\ ■ K ) T\ + (1 «- 2) . (4) 
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The exponential term represents the plane wave of the pion field. The structure of the wAN 



coupling follows from the usual interaction lagrangian [13 in the non-relativistic reduction. 



(T*) is the spin (isospin) transition matrix [XT] for the A excitation, and r = f\ — r*2 is the 
relative coordinate of the two nucleons in the deuteron. The symbol (1 •*->• 2) in eq. (|J) shall 
indicate that the operator has to be properly symmetrized with respect to the nucleonic 
coordinates, i.e. each of the two nucleons in the deuteron may be excited to a A resonance. 



Following refs. [TgJTjJ], the quantum mechanical state 

Ipa;v) = Fify m (5) 

is called the uncorrelated source function for the AN intermediate system, with \%pd) being 
the deuteron wave function. After propagation of the AiV system and absorption of a second 
pion, one obtains the source function for the AA state, 

\p aa) = Fl(k b ) G AN Fl(k a ) \tfj d ). (6) 

In our model, the full propagator Gan is given by 

AN can + |Ta(sa) — Tan — Van ^ ^ 

It contains the excitation energy €an, the energy-dependent width Ia(sa) °f the resonance 
T9fl , the operator of the kinetic energy Tan and the interaction potential Van of the inter- 



mediate AiV system. €an and sa are fixed by imposing energy conservation at the 7riVA 
vertex. The potential Van is constructed in a meson exchange model with tt, p, uj and o 
exchange taken into account. A diagrammatic representation of the residual AiV interaction 
is given in fig. |2| (a,b). For more details, we refer the reader to refs. [|i4| , |T3| where also explicit 
expressions for Van can be found. 

In eq. (|B|) we have not yet considered that the AA source function has to be symmetric 
under exchange of the pions a and b. According to the value T of total isospin in the im 
system, we define the symmetrized operator 

Fl(k a , h) = ^ {Fl{h) Gan F^fy ± F^fy Gan (8) 



, + if \TM T ) = |00) = -4=|7T + 7T- - 7T 7T + 7T 7T + ) 

with <( V3 

- if |TM T ) = |10) = ^|vr+7T- - 7T-7T+) 

The complete matrix element for the AA mechanism can now be written as 

Maa = (ip n \ (4> P \ V A a~*nn Gaa Fl w (k a , k b ) \ip d ) (9) 

Here, ip n ( p ) denote the distorted waves of the projectile neutron (target proton) which are 
calculated in eikonal approximation. In analogy to eq. (|7|), the full propagator for the A A 
system is given by 

e AA + | [r Aa (sAj + r A> (sAj] - T AA - Vaa ' ^ 

Especially, this propagator includes the residual interaction Vaa m the excited AA system, 
see also fig. 0(c). Vaa, as well as the transition potential Vaa->nn between the intermediate 
AA system and the NN state, are both constructed within the same meson exchange model 
as Van- To the A A — > NN transition potential, of course, only the tt and p meson exchanges 
contribute. The AA excitation can thus be separated into a spin-longitudinal part (with 
7r-like coupling to NN) and a spin-transverse part (with p-like coupling to NN). As we 
will demonstrate later, this spin-structure of the A A excitation is important for the angular 
distribution of the pions. 

The matrix element corresponding to the second mechanism where only one A is excited, 
fig. ||(b), is found to be 

Man = (Vnl FL$a, h) \ip d ) (11) 

Here, the operator F^ n is obtained from by simply replacing for the second pion and 
with their hermitian adj ungates S and T, respectively. 

In order to obtain the 27r production amplitude, the coherent sum of the two contributions 
of eqs. (|9|, [TT]) has to be taken. Due to the fact that the residual interactions Van and Vaa 
couple states of different quantum numbers, the concrete evaluation of the matrix elements 
leads to a system of coupled integro-differential equations. It turns out that they can be 



solved in a very effective way with the so-called Lanczos method | 20| . This procedure has 
already been described in ref. |H| and can easily be applied to the present calculation. 



III. RESULTS AND DISCUSSION 



A. Parameters of the model 



Input parameters of the model are the meson and baryon masses, coupling constants, 
and formfactor cutoffs at each vertex. The values we used in our calculations are given 
in table ||[ Most of these quantities, of course, are physical observables that are already 
determined from other experiments, e.g. from NN scattering [^l]]. Only the sigma meson 
and cutoff parameters remain to be fixed. They have been adjusted to fit the experimental 
data for the one pion production N + N —>- 2 }Itt. We made sure that not only the total cross 
section but also the angular distributions and analyzing powers are correctly reproduced in 



the A resonance energy region [22]. After this consistency check, we may now continue and 



discuss our results for the 27r production n + p ^ 2 H (jrir). 



B. The deuteron spectrum 

In fig. |3|, we present the experimental deuteron spectrum d 2 a / ' dk^dVL^ as it is measured in 
the reaction n + p — > 2 H(tttt) at a beam momentum k n = 1.88 GeV and in forward direction 
6d = 0°. The production of both neutral (7r°7r°) and charged (7r + 7r~) pion pairs is possible. 
The available phase space for the 2n production is indicated by the dashed line. 

Obviously, the cross section does not follow the phase space but shows a characteristic 
structure which three pronounced peaks. The two outer maxima correspond to a kinematical 
situation where the invariant mass M n7T of the 2n system is minimal, i.e. where = Im^. 
The broad central maximum, on the other hand, is located at the position where M 7Tn has 
its maximal value (511 MeV in the present case). In the CMS, this situation will be realized 
if all the kinetic energy is taken by the two pions and the deuteron remains at rest. 
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As can be seen from fig. |], our model is able to explain the characteristic structure of 
the experimental spectrum. The solid line represents the full result with all intermediate 
interactions included. The A A mechanism, for which the Feynman diagram was presented in 
fig. [j](a), is in fact the only important contribution to 2ir production in the energy range under 
consideration. The cross section contribution due to the A7V" mechanism of fig. [TJ(b) is almost 
negligible. Since the optimal energy for this second mechanism would be y/s = MA+M—m n , 
its suppression in the ABC energy region is not surprising. 



With focus on the A A excitation and following the arguments of ref . [JLT] , the origin of 
the different maxima can be understood quite easily. Let us define 

K = p a + Pb and k = p a - p b , (12) 

where p a , pb denote the four-momenta of the two pions. The excitation of a AA system 
is most effective if the invariant masses of the two A's are approximately equal (and hence 
both close to the resonant mass). In this case we have 

s Aa = \ (Pd + K + k) 2 = 1 (p d + K - k) 2 = s Ab , (13) 

and since K • k = this is equivalent to 

Pd-k = 0. (14) 

This condition can be fullfilled in two ways: 

1 . k = and therefore M W7T = VK 2 = 2m- K , or 

2. f>d = in the CMS of the two pions (which means that the deuteron restframe and the 
2tt restframe are identical). In the latter case, the two pions carry the whole kinetic 
energy, and M nn = max. 

The first situation leads to the outer maxima and corresponds to the parallel decay of 
the two A's because of p a ~ pb, i.e. the relative momentum is very small. The second 
situation explains the central maximum and corresponds to the antiparallel decay of the 
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AA excitation for which the relative momentum is maximal. In the deuteron restframe, the 
different kinematical configurations for the three maxima can be visualized as depicted in 
fig- | 



C. Angular distribution of the pions 

The kinematical situations leading to the three maxima of the deuteron spectrum are 
also clearly visible in the angular distribution of the pions. In fig. |5], we show the triple 
differential cross section d 3 a/dk ( idQ ( idQ 7T for different deuteron laboratory momenta. It is 
plotted as a function of n which is the angle between the relative momentum k of the two 
pions and the beam axis, as given in the deuteron rest frame (see also fig. [5]). 

For the two outer maxima (k l ^ h = 1.1 GeV and k x ^ h = 1.9 GeV), the pion momenta are 
nearly parallel and hence the relative momentum is dominantly perpendicular to the beam 
axis. Consequently, the differential cross section reaches its largest value at cos 6^ = 0. 

For the central maximum (/cjf b = 1.5 GeV), however, the pion momenta are antiparallel, 
and all angles 8^ are kinematically possible. If there was no spin dependence of the excitation 
and the residual interaction, an isotropic distribution would be the result. The observed 
angular variation of the cross section is thus reflecting the spin structure of the process 
which mainly follows from the 7r and p exchange contributions to the NN — > A A transition 
potential. 

In order to examine this spin structure in more detail, we will neglect for the moment the 
residual interaction and the possibility of spin-flips. Then, the spin-longitudinal 7r exchange 
leads to an operator proportional to 

'(Si ■ q) (S\ ■ k a )] [(S 2 ■ q) (S\ ■ k b )] ~ (q ■ k a ) (q • k b ) , (15) 

where k a , k b are the pion momenta and q is the momentum transfer from the neutron pro- 
jectile to the proton target. If both the pion momenta are parallel (or antiparallel) to q, the 
corresponding cross section will be largest. This results in a maximum at cos 6^ = 1 in the 
spin-longitudinal channel. 



In complete analogy, the spin-transverse p exchange has the operator structure 

"(£ x q) (S\ ■ k a ] [(S 2 x q) (S\ ■ k b )] ~ (qx k a ) (qx k b ) . (16) 

Therefore, the maximum in the spin-transverse channel is reached if both the pion momenta 
are perpendicular to q, i.e. if cos 9 n = 0. 

Indeed, the theoretical curves in fig. [7| exactly reflect the expected shapes. We conclude 
that a measurement of the pion angular distribution would be very helpful in order to reveal 
the spin structure of the interaction. 

D. Influence of the residual interaction 

Effects of residual interactions can be examined in both the AA and the AiV system since 
these are the intermediate configurations in our model. The influence of the corresponding 
interaction potentials on the theoretical spectrum is shown in fig. |8| for the case of forward 
scattering (9d = 0°). 

The dash-dotted line was calculated without any residual interaction, i.e. with V = 
for the intermediate states. The typical ABC structure with its three peaks as discussed 



in section [II IB] is clearly present. Our result comes quite close to the fully relativistic 
calculation of Bar-Nir et al. [§ which also does not include the residual interactions. 

By choosing V = Van but still neglecting the direct AA interaction, we obtain the 
dashed line. Here, the maxima are even more pronounced. The enlargement of the peak 
cross sections is due to the reduction of the excitation energy for the AiV system which is 
caused by the attractive AiV potential. For y/s < 2M& as in the present case, the whole 
system is thus closer to the resonance energy. The relative position of the two invariant A 
masses remains however unchanged since in our model they are fixed independently of Van- 

On the other hand, the AA potential can directly influence the relative wave function 
of the AA system via a redistribution of momentum within the two particle system. The 
solid line in fig. R] has been calculated with the full residual interaction V = Van + Vaa- 
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The inclusion of Vaa results in a transfer of strength from the maxima to the regions 
between which are kinematically less favoured. This transfer takes place because the optimal 
configuration of equal A masses can now be reached even if the initial energy distribution 
of the two pions was asymmetric. As a consequence, the cross section is enhanced between 
the peaks and simultanously the magnitude of the peaks is reduced. This reduction is more 
prominent for the outer maxima since they are related to a small relative momentum of the 
AA system. For the central maximum, the relative momentum is large and therefore only 
a slight influence of the AA potential is observed. 

Obviously, the residual interactions play quite an important role in the 2tt production 
and significantly influence the deuteron spectrum. The inclusion of these effects in the 
calculation clearly improves the agreement between theory and experiment. 

E. Angular and energy dependence of the deuteron spectrum and the total cross 

section 

Studying the dependence of the deuteron spectrum on the scattering angle 9d, we find 
that the theoretical cross section for A A excitation is decreasing too fast. Exemplatory 
we present in fig. |9] our calculations for 9d = 4.5° and 9d = 7.5°. The experimental data 
0] are underestimated by a factor 2 ~ 3. Other microscopic calculations, e.g. assuming 
a two nucleon exchange as the dominant reaction mechanism, yield a comparable angular 
dependence and are also not able to describe the data [PJlCfl. The reason for this behaviour 
is not understood. Maybe more than only one production mechanism has to be taken into 
account in order to solve this problem. 

Fig. ^ is demonstrating the energy dependence of the deuteron spectrum in our model. 
One recognizes that the central maximum is getting more pronounced if the total energy 
increases. We remind the reader that this peak corresponds to the AA excitation with 
maximal relative momentum, which already gives the natural explanation of the observed 
effect. 
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After integration over deuteron momentum and scattering angle, the total cross section 
for the n + p — > 2 H(irir) reaction is obtained. In fig. |11|, our theoretical result is compared 
to the experimental data In the A A resonance energy region, i.e. for a neutron 

momentum of k n w 2 GeV, the agreement is quite good. This confirms our assumption 
that the 2tt production in the regime of the ABC effect is dominated by the AA excitation 
mechanism. On the other hand, one recognizes that the experimental cross section close 
to the 2-7T production threshold (k n = 1.19 GeV) is understimated by nearly two orders 
of magnitude. In this case, the total energy of \fs = Md + 2m n <C 2Ma is simply not 
sufficient for the excitation of a AA system, and other production mechanisms will be 
more important. To mention just two possibilities, the pions could be produced by s-wave 
rescattering or via excitation of the A* (1440) resonance, as recently discussed in p5|| . Since 
there is also experimental evidence p4| that the ABC structure completely disappears when 
approaching the 2n production threshold, further theoretical investigations of the low energy 
region would be of high interest. 



IV. SUMMARY AND CONCLUSIONS 

To summarize we have shown that the AA excitation is the dominant reaction mechanism 
in the n+p — > 2 H (tttc) two pion production at k n = 1.88 GeV/c. It is able to explain the ABC 
effect as observed in the experimental deuteron spectrum. Hereby, the residual interactions 
in the intermediate AA and AA states play an important role. We obtain good agreement 
with experimental data at forward scattering but too fast a decrease of the cross section 
for higher scattering angles of the deuteron. The total cross section is also underestimated 
close to the 2tt production threshold where other reaction mechanisms get more important. 
Furthermore, we found that the spin structure of the AA excitation directly influences the 
angular distribution of the two pions. We conclude that the 2tt production in the ABC 
energy region may well serve as a tool for closer examination of direct AA interaction and 
transition potentials. 
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TABLES 



TABLE I, Parameters of the model. 

Unn/^ gWjTT /^aa/4vt A a [GeV] m a [MeV] 

7T 0.081 0.32 0.0031 1.1 138 

p 5.4 21.6 0.286 1.4 770 

w 8.1 a 8.1 a 1.7 783 

a 6.9 a 6.9 a 1.6 580 

a 5aiV7v/ 47r res P- 5aAA/ 47r is g iven - 
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FIG. 1. Feynman diagrams of 2n production via A excitations in the n+p -^ 2 H (tttt) reaction: 
(a) AA excitation, (b) AiV excitation. The shaded areas symbolize the residual interaction in the 
intermediate AA and AiV states, respectively. 
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FIG. 2. Residual interactions in the meson exchange picture, (a) Direct and (b) exchange 
contribution of the A./V potential, (c) AA potential. The mesons taken into account are ir, p, to, 
and a. 
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FIG. 3. Experimental deuteron spectrum d 2 a / 'dk^dfl^ of the reaction n+p -^ 2 H (tttt), measured 
at k n = 1.88 GeV and = 0°. Data have been taken from Qj. The dashed line is the phase space 
for the 2n production (in arbitrary units), and M wn is the invariant mass of the 2tt system. 
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FIG. 4. Theoretical result in comparison with the experimental spectrum of fig. ||. Solid line: 
full calculation with AA and AiV mechanism included; dashed line: only AN excitations, (see 
also fig. H). 
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FIG. 5. Kinematical configurations in the deuteron rest frame pa = 0. Right and left: Momen- 
tum configuration for the outer maxima of the deuteron spectrum; middle: momentum configura- 
tion for the central maximum. 
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FIG. 6. Angular distribution ^a/dkddtlddQn of the pions, for k n = 1.88 GeV and 6d = 0°. 0^ 
is the angle between the relative pion momentum and the beam axis in the deuteron rest frame. 
The different deuteron laboratory momenta of 1.1, 1.5 and 1.9 GeV correspond to the left, central, 
and right maximum of the deuteron spectrum. 
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FIG. 7. Contributions of the spin-longitudinal (LO) and spin-transverse (TR) component of 
the transition potential NN — ► A A to the differential cross section at kd = 1.5 GeV (without 
residual interaction). 
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FIG. 8. Influence of the residual interaction on the deuteron spectrum, for k n = 1.88 GeV and 
9d = 0°. Solid line: full calculation with both AiV and A A potential; dashed: with AiV potential 
but without AA potential; dash-dotted: without any residual interaction. 
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7.5°, respectively. The solid line is the theoretical result. Experimental data from ref. H. 





3.0 


i 1 i 1 i 1 




0.30 


i 

W 


2.5 


- k n =2.5 GeV 




0.25 


> 






> 




CD 


2.0 




CD 


0.20 




1.5 














0.15 


"D 

a 


1.0 




"D 

a 

T3 


0.10 


~D 

c\j 


0.5 




T3 

CM 


0.05 


"D 


0.0 




T3 


0.00 



0.8 1.2 1.8 2.2 2.8 
k d [GeV] 



k =1.46 GeV 



0.8 




k d [GeV] 



1.6 



FIG. 10. Theoretical deuteron spectra in our model for forward scattering 0^ = 0° at beam 
momenta k n = 2.5 GeV and k n = 1.46 GeV, respectively. 
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FIG. 11. Total cross section for n + p -^ 2 H(7T7t) as a function of the neutron momentum. 



Experimental data from j|||,^2.24]. For k n = 2.09 GeV, we have = 2Ma- The threshold for 
2tt production is at k n = 1.19 GeV. 
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